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Abstract 

The neutron charge radius is studied within a light-front model with different spin coupling 
schemes and wave functions. The cancellation of the contributions from the Foldy term and Dirac 
form factor to the neutron charge form factor is verified for large nucleon sizes and it is independent 
of the detailed form of quark spin coupling and wave function. For the physical nucleon our results 
for the contribution of the Dirac form factor to the neutron radius are insensitive to the form of 
the wave function while they strongly depend on the the quark spin coupling scheme. 



PACS numbers: 12.39.-x,13.40.-f,13.40.-Gp,14.20.-c 



I. INTRODUCTION 



Nowadays, there is a renewed interest in the nucleon electromagnetic form factors due to 
recent precise experiment sQ]. The electroweak static observables are receiving theoretical 
attention as well, in particular, the neutron charge radius was recently studied with respect 
to its relativistic origin . From the nonrelativistic point of view the nucleon wave 

function with SU(6) symmetry implies a zero neutron charge radius and \i v j \i n = —m^j (2m) 
(mjy and m are the masses of the nucleon and constituent quarks, respectively). However, 
the quark spin-spin interaction that implies the nucleon-delta mass splitting allows the 
dynamical breaking of the SU(6) symmetry and a nonvanishing neutron charge mean square 
radius (r 2 ). The flavor identical quarks are pushed out, while the u quark stays around the 
neutron center, which gives a negative charge mean square radius, in agreement with the 
observed sign p. 

The experimental value r 2 n = —0.113 ± 0.005 fm 2 J(| is coincidentally near the contri- 
bution of the Foldy term (f^f- = -0.126 fm 2 ), and for that reason the contribution of the 
Dirac form factor (Fi n ) is quite small, r\ n = 0.013 ± 0.005 fm 2 . The electric and magnetic 
form factors (Sachs form factors) are given by: 

G EN {q 2 ) = F lN {q 2 ) + -^ T F 2N {q 2 ) , 

4777, 

G MN (q 2 ) = F 1N (q 2 ) + F 2N (q 2 ) , (1) 

where N = norp, F 2 n is the Pauli form factor, q^ is the momentum transfer. The magnetic 
moment is fi^ = Gmn(0) and the charge mean square radius is r 2 N = 6 dGE d p q ^ \ q 2 =0 . 

The naive physical picture of the spin-spin interaction dominating the negative neutron 
charge radius is confronted with the fact that the recoil effect from the Foldy term dominates 
the charge radius. Thus, if one associates the intrinsic charge distribution with F ln , one 
comes to the wrong conclusion that the spin-spin interaction is not relevant for r 2 . The 
recent work of Isgur[2] clarified this issue, the rest frame charge distribution of the neutron 
should be associated with the form factor Gg„ and not with F ln . He found that in the the 
first relativistic correction to the Dirac form factor cancels the Foldy term in a relativistic 
model of the nucleon within a quark-diquark picture. Consequently, the nonzero value of 
the neutron charge radius should reflect the intrinsic charge distribution and the detailed 



quark spin dynamics. In fact, in Ref. 
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it was shown that the different forms of relativistic 
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quark spin coupling in the nucleon have dramatic effects on the neutron charge mean square 
radius. The relativistic spin coupling coefficients, which depend on momentum, effectively 
lead to the breaking of the SU(6) symmetry as discussed in Ref. ^], and to a nonzero neutron 
charge form factor. 

In this work, we focus our attention on the neutron charge mean square radius and our 
aim is to explore, within a light-front relativistic model| , the detailed dependence of rf n , 
obtained from iq n , on the nucleon size that is parametrized by the proton charge radius 
(r p ). The proton charge radius controls the relativistic or nonrelativistic nature of the 
constituent three-quark system. For r p ^> m , the inverse constituent quark mass, the 
nucleon approaches a nonrelativistic system of quarks. On the other hand for r p ~ mT l ~ 
1 fm (the real nucleon size and constituent mass scales), the three-quark system demands a 
relativistic description of the internal quark motion. The parameters of the light-front model 
are adjusted to different nucleon sizes, with several forms of momentum component of the 
wave function and quark spin coupling schemes. By changing the parameters we are able 
to shift smoothly from relativistic to nonrelativistic regimes. In the nonrelativistic regime 
we show quite generally the cancellation between r\ n and the Foldy term contribution to 
the neutron charge radius, while at the physical nucleon scale the value of r\ n is strongly 
dependent on the choice of quark spin coupling scheme. 

Following our previous work[3j, where we have studied the nucleon electroweak form 
factors using different forms of relativistic spin coupling between the constituent quarks to 
form the nucleon, we use an effective Lagrangian to describe the coupling of the quark spin. 
The form factor calculation keeps close contact with covariant field theory. The starting 
point is the impulse approximation for the nucleon virtual photon absorption amplitude, 
which is projected on the three-dimensional null-plane hypersurface, x + = x° + x 3 = 0, (see, 
e.g., Ref. |7(). The three-dimensional reduction is done by integrating over the individual 
light-front energies (k~ = k° — k 3 ) in the two-loop momentum integrations of the impulse 
approximation. The relative light-front time between the particles is eliminated in favor of 
the global time propagation [8j. Then, the momentum component of the nucleon light-front 
wave function is introduced into the remaining form of the two-loop momentum three- 
dimensional integrations which define the matrix elements of the electroweak current 
llo| . In general, the vertex function depends on the quark momentum variables Xi, and is 
chosen to be totally symmetric. We make the common assumption that the vertex function 
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depend on the free three-quark light-cone energy, Mq (Eq. 9), which is the simplest totally 



£ N _ 3q = am N e lmn ^ il) ir 2l ^f m) ^ in) ^ N + (1 - a)e lmn ^ {l) ir 2l ^f m) ^ in) zd^ N + h.c. (2) 



where t 2 is the isospin matrix, the color indices are {I, m, n} and e lmn is the totally anti- 
symmetric symbol. The conjugate quark field is = C^ T , where C = i'j 2 ^ is the charge 
conjugation matrix; a is a parameter to choose the spin coupling parameterization. 

In Ref. P| we have tested different spin couplings for the nucleon in a calculation of nucleon 
electroweak form factors. We have found that the neutron charge form factor constrains the 
relativistic quark spin coupling schemes. The neutron data below momentum transfers of 
1 GeV/c suggested that the scalar pair (a = 1) in the effective Lagrangian is preferred. 
In that study, Gaussian and power law momentum components of the wave functions were 
used. Here we enlarge the set of momentum components of the wave functions to allow a 
wide variation of parameters, to shift from relativistic to nonrelativistic regimes. 

This work is organized as follows. In section II, it is given a brief description of the 
macroscopic and microscopic forms of the nucleon electromagnetic current appropriate for 
the light-front calculations. In section III, we present the numerical analysis of the static 
nucleon observables for different model assumptions. A summary and conclusion is presented 
in section IV. 

II. NUCLEON ELECTROMAGNETIC CURRENT 
A. Macroscopic matrix elements 

The eletromagnetic form factors are extracted from the plus component of the current 
for momentum transfers satisfying the Drell-Yan condition q + = q° + q 3 = 0. The contribu- 
tion of the Z-diagram is minimized in a Drell-Yan reference frame while the wave function 
contribution to the current is maximized mum 

. In particular, the Breit-frame is 
chosen, with four momentum transfer q = (0, <fj_,0), such that (q + = 0) and q± = (g 1 ,? 2 ). 
The nucleon momentum in the initial state is p = (v/4 + rn 2 N , —^-,0) and in the final state 



symmetric scalar function of the quark momentum variables. 
The effective Lagrangian for the N — q coupling is written as 



a, 



is given by p' = (J 2j- + ^ u_ 
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The macroscopic matrix element of the nucleoli electromagnetic current J^(q 2 ) in the 
Breit-frame and in the light-front spinor basis is given by: 

(sVUq 2 )\s) = u(p',s') (iW) 7 + + l a -^F 2N {q 2 )^u{P,s) 

= ^-(s'\F 1N (q 2 )+^^^<l± ■ {<? x n)\s) , (3) 

where and F 2 n are the Dirac and Pauli form factors, respectively, while n is the unit 
vector along the z-direction. 

The light-front spinors are defined as: 

/ 



2y/p + m N 



v Pauli 



\ 







(4) 



and the Dirac spinor of the instant form is given by 



p' + rriN 

u D {p, s) -- 



v Pauli 



V 



o 



(5) 



'2m(p° + m 

which carries the subscript D. The Melosh rotation is the unitary transformation between 
the light-front and instant form spinors, which is given by: 

[Rm( P )U = (*1 P+ * m ~'*f 7\ > = *»(P. *) • (6) 

B. Microscopic matrix elements 

The microscopic matrix elements of the nucleon electromagnetic current is derived from 



the effective Lagrangian, Eq.(J2J), within the lig 
resented by the two-loop diagrams of figure 



it-front impulse approximation which is rep- 
The complete antisymmetrization of the 



n 



quark states implies that the matrix element of the current is composed by four topologically 
distinct diagrams depicted in the figure. The matrix elements of the electromagnetic current 
are calculated considering only the process on quark 3, due to the symmetrization of the 
microscopic matrix element after the factorization of the color degree of freedom. The four 
distinct current operators Jt^, j3 = a,b, c, d, are constructed from the Feynman diagrams of 
figure la to Id, respectively |3J]. 

The microscopic operator of the nucleon electromagnetic current, Jjy, is the sum of each 
amplitude represented by the diagrams (la) to (Id): 

= JM) + 4<4W) + 2 Jt N {q 2 ) + 2J+ N (g 2 ) ; (7) 



where the weigthing factors comes from the identity of quarks 1 and 2, and another factor 2 
multiplying J^ N comes from the exchange of the pairs in the initial and final nucleons, which 
gives the same matrix element as a consequence of time reversal and parity transformation 
properties. 

The light-front momentum are defined as k + = k° + k 3 , k~ = k° — k 3 , k± = (k 1 , k 2 ). 
In each term of the nucleon plus component of the current, from J^ N to j£ N , the quark 
momenta are on-fc~-shell. The total plus and transverse momentum components of the 
intermediate states satisfy conservation laws. Thus, the components of the momentum kf 
and k~2 are bounded, such that < k~l < p + and < k 2 < p + — kf jl^. 

The two-loop Feynman diagram of figure la corresponds to 

(sUW)i*> = 2 P +2 (N\Q q \N) f ^ 1± dkt^dl4 0{p+ _ kt)e{p+ _ ^ _ K) 

J fv^ ^3 

Tr [(^ 2 + rn) (am N + (1 - a)f) (fii + m) (am N + (1 - a)fi)] 

u(p', + m))7 + (^ 3 + m)u(p, s)tt (M' 2 )* (M 2 ) , (8) 

where k\ = m 2 and k\ = m 2 . The momentum component of the wave function is \I/(Mq) 
and the free three-quark squared mass is defined by: 

Ml = + ^-^! + - P l , (9) 

Ki ft 2 % 

and Mq 2 = Mq (k 3 — > k f 3 , pj_ — > The electromagnetic quark current operator is ^Qg^^f, 
with Q 9 the charge operator and \I> the quark field. 

The other terms of the nucleon current, represented in figures lb to Id, are written below: 

d 2 ki±dkid 2 k 2 _Ldk2 



(AJtstf)\*) = P +2 (N\Q q \N) f 1 2 ^ + - k ^ P+ - k l - k2 

J nJ-j^ rC<2 nJg 

u(p', s')(jk' 3 + m)7 + (^ 3 + m) (am N + (1 - a)f) (^i + m) 



x (am N + (1 - a)p") (fe + m)u(p, ^(M^^M 2 ) , (10) 
(sV+ N (d 2 )\s) = P +2 (N\r 2 Q q r 2 \N) [ * k ^fJ^ d ® 0(p+ - k+)6{p + - # - H) 

J \ 2 3 

u(p', s')(fti + m) (am N + (1 - a)f) (jt s + m)^ + (f,' 3 + m) 

x (am N + (1 - a)p") + m)w(p, s)#(Mo 2 )#(M 2 ) , (11) 

J rCi /Co rCq 



Tr 



Jc+U+U- 
h l K 2 ^3 

(am N + (1 - a)f) (jt' 3 + m)7 + (^ 3 + m) (am N + (1 - a)fy (^ + m) 



s')(fe + m)«(p, s)^(M' 2 )^(Ml) . (12) 



In our study, for simplicity and as is explicit in Eqs.(JHJ to (|12|h the same momentum 
wave function is chosen for both N — q couplings. The momentum part of the wave function 
in the microscopic matrix element of the current is chosen from different models as we will 
show in the next section. 



III. RESULTS FOR STATIC NEUTRON OBSERVABLES 

In this section we present our theoretical results for static nucleon charge square radius 
assuming the dominance of the lowest light-front Fock state component in the nucleon wave 
function corresponding to three constituent quarks. By itself this is a strong constraint on 
the static observables and essentially the results are mostly dependent on the constituent 
quark mass and scale with proton charge radius, as it will be shown. We use a constituent 
quark mass value of m = 0.22 GeV from Refs. P. . 

Within the above assumptions, we show the effects of different relativistic spin couplings 
and momentum wave functions of constituent quarks for the neutron and proton charge 
square radius. The key point in our discussion is the detailed dependence of 

on the nucleon size, parameterized by the proton charge radius (r p ), which controls the 

relativistic or nonrelativistic nature of the constituent three-quark system. The Foldy term 

contribution to the neutron charge square radius is 

2 3 /i n 

r2n ~ 2^ ' 

and the square of the neutron charge radius is the sum of both contributions, i.e., r\ = 

y 2 I 71 2 

' In. ' 'In- 

The parameters of the light-front model are changed to allow different nucleon sizes for 
different forms of momentum component of the wave function and quark spin coupling 
schemes. By modifying the parameters a smooth shift from relativistic to nonrelativistic 
regimes is obtained. The nonrelativistic regime of the constituent quark system is charac- 
terized by r 2 > 1 fm 2 (about the inverse of the constituent quark mass squared), and the 
relativistic regime by r 2 < 1 fm 2 (the real nucleon size and constituent mass scales). 

The correlations between the neutron charge square radius with magnetic moment and 
proton charge square radius are investigated with a different momentum part of the nucleon 
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light-front wave function for each quark spin coupling scheme. Among the observables, the 
neutron charge radius plays a special role; its correlation with the magnetic moment depends 
on the quark spin coupling scheme^. 

A. Model wave functions 

The matrix elements of the microscopic nucleon plus component of the current, Eqs.(jHJ) to 
(|12)1. are evaluated for different A — q spin couplings and momentum part of the wave func- 
tion. The different models of the momentum component of the wave function corresponds 
to the choices of the harmonic and power-law forms 
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14, 



*ho = A HO exp(-M 2 /2/3 2 ) , ^ Power = A Power (l + M*/{3 2 y P (13) 



and modified harmonic and power-law wave functions 



exp(-M 2 /2/3 2 ) (1 + MHPyv 

fll — Mq ' Power ~~ Power (31 — Mq 



The normalization is determined by the proton charge. The width parameter is (3. The 
free mass of the three-quark system satisfies Mo > 3m and (3\ has to satisfy the constraint 
Pi < 3m to avoid scattering poles in the bound state wave function of Eq.(jHJ). This property 
is consistent with color confinement which prevents scattering states of three quarks to be 
relevant. 

The wave function models of Eq. ()14j) are inspired in the general form of the lowest 
Fock-state component of the nucleon wave function in QCD light-front field theor y w here 
the complete wave function is an eigenstate of the complete mass operator equationjlll. llfil ] . 
The lowest Fock state component of the hadron wave function satisfies an effective mass 
operator equation for constituent quark degrees of freedom, in which the effective interaction 



contains in principle all the complexity of QCD 
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The general form of lowest Fock 



component in terms of the constituent quark degrees of freedom has the term {(3\ — Mq) -1 , 
where (3\ plays the role of the mass of a bound system. Here, we use the models for the 
wave function from Eq. (pPfj) just to enlarge our possible choices of momentum components 
of the wave function, while still keeping connection with the basic QCD theory. 



S 



From general QCD perturbative arguments a power-law fall-off with p = 3.5 is predicted 



for ^Power 







14 1 . The correlations between static electroweak observables are not sensitive 
to p as long as p > 2 and we choose for our calculations p = 3 in both \l/p OW e r and 

^Power- For large virtualities \l/power ~ M ~ 2p and ^'p owcr . ~ MJ~ 2 " ,+1 \ with p = 3, which is 
above and below the QCD power fall-off, respectively. As we are going to show below these 
different assumptions for the wave function have only little effect in the correlations between 
the static observables. 



B. Numerical Results 

In Figs. 2 to 5 we show results for the correlations between static neutron charge square 
radius, magnetic moment and proton radius. Our calculations are done for different spin 
couplings of quarks, i.e. a = 0, 1/2, 1 in the effective Lagrangian of Eq.(J2J), and momentum 
wave functions of a harmonic oscillator (HO) (Gaussian) and a power-law (Power) form 
(p = 3) from Eq. (fT3*j) and the modified forms from Eq. (|14|1. 

The correlation of the static observables is found by variation of the (3 and (3\ parameters. 
In the Gaussian wave function of Eq. (|T3*j) two limits are noteworthy, f3 — ► leads to an infinite 
size of the nucleon corresponding to the nonrelativistic limit and (3 — > oo is the zero radius 
limit corresponding to the strong relativistic limit. In the power-law wave functions for 
(3 — > oo the relativistic limit is approached. However in this case of (3 — > one does not 
approach the nonrelativistic limit because the typical momentum scale of the wave function 
is the quark mass and not zero. The modified Gaussian and power-law wave function of Eq. 
()14j) are much more flexible, because the nonrelativistic limit can be approached by taking 
f3\ — > 3m, which corresponds to zero binding energy, and the quark system swells to infinity. 
Taking (3 going to infinity in the Gaussian and in the modified Gaussian models, the nucleon 
tends to zero size, and to the extreme relativistic regime. In our results we explore a wide 
range of values of (3 and j3\. 

In Figure 2 results are shown for the neutron charge radius as a function of the neutron 
magnetic moment for a = 0, 1/2, and 1 as well as HO, Power, modified HO and Power 
momentum wave functions. The results are quite insensitive to the different shapes of 
the momentum wave functions, however strongly dependent on the quark spin coupling 
as we have already found in Ref. j^j]. The present extension to several different forms of 
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wave function confirms the previous findings, and moreover the confining or nonconfining 
behavior of the wave function is not important for the neutron radius as long as the magnetic 
moment is fitted. The nonconfining feature of the wave function does not change our previous 
conclusion, i.e., the experimental data for the neutron charge radius favors the scalar coupling 
between the quark-pair, while the gradient spin coupling (a = 0) completely disagrees with 
the experimental value. 

As our main point is to study the neutron radius as a function of the nucleon size, 
parameterized by the proton charge radius, next we show in Figure 3 the neutron charge 
square radius as a function of the proton charge square radius, for the same set of calculations 
presented in Figure 2. The different models of quark spin couplings (for a equal to 0, 1/2 
and 1) are shown in Figure 3 and represent a systematic behavior that once again is quite 
independent of the form of the momentum component of the wave 



unction. For the chosen 

constituent mass (m = 0.22 GeV) the experimental points from [17|, 18, 3] data are within 
the width for the theoretical results for the scalar coupling. The results for the modified 
Gaussian and power-law wave functions are obtained by changing either (3 or f3\ in Eq. (fT^|) . 
which leads to a small spread of the results seen in the figure. 

The functional dependence of the individual contributions r\ n and r\ n to the neutron 
charge square radius with proton charge square radius is shown in Figure 4. The relativistic 
and nonrelativistic regimes are identified in the figure. For the proton charge square radius 
below 3 fm 2 , we observe the intrinsic relativistic behavior of the quark motion through the 
wide separation of the results for r\ n obtained for different quark spin coupling schemes. 
The nonrelativistic regime is seen for r 2 above 3 fm 2 where the calculations tend to be flat 
and with magnitude which cancels to some extent the contribution of Foldy term, r\ n . At 
the physical nucleon size scale r p ~ 0.8 fm, the quark motion is quite relativistic, where r\ n 
has a sensitive dependence on the spin couplings, i.e. a. The result for r\ n is about zero 
for the scalar coupling which implies in this case the dominance of the Foldy term in the 
neutron charge radius. The dependence on the details of the momentum component of the 
wave function is quite small. 

To be complete, we present in Figure 5 the function defined by r 2 (r 2 ) over the same range 
of proton sizes as in Figure 4. The results for the neutron square charge radius show a smooth 
trend in all models and spin coupling schemes from the physical scale (relativistic regime) 
towards the nonrelativistic regime. The cancellation of the Foldy term in the nonrelativistic 
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regime leads to the small values of r\ for all models investigated despite the differences in 
the quark spin coupling schemes and wave functions. We show the cancellation between 
r\ n and the Foldy term contribution, r\ n to the neutron charge radius in the nonrelativistic 
regime for the particular coupling of Eq. (J2J), while at the physical nucleon scale the value 
of r\ n is strongly dependent on the choice of quark spin coupling scheme. 

IV. SUMMARY AND CONCLUSION 

In this work we have studied in detail the effect of the nucleon size scale in the individual 
contribution of the Dirac form factor to the neutron charge square radius, using a rela- 
tivistic light-front model with constituent quarks. The model is constructed with different 
relativistic spin coupling schemes and wave functions. The wave function parameters were 
adjusted to different nucleon sizes, parameterized by the proton charge radius, which allowed 
us to investigate the cancellation between the contributions from the Foldy term (r| n ) and 
Dirac form factor (rf n ) in the neutron charge square radius that occurs in nonrelativistic 
regimes Q]. 

First we extend the previous analysis of the neutron static observables using more 
general forms of the momentum component of the wave function, which have a nonconfining 
tail and we find a model independent correlation of the Foldy term and r\ n with the proton 
charge radius for each spin coupling scheme. The existence of the model independent cor- 
relation of r\ n and r\ n with the proton charge radius allowed the study of the nucleon size 
scale effects on a physical ground. 

Our calculations show that the cancellation between r\ n and the Foldy term contribution, 
r\ n to the neutron charge radius happens indeed to a large extent in the nonrelativistic 
regime, independent of the detailed form of quark spin coupling scheme and wave function. 
The nonrelativistic regime is seen in our calculations for the square of the proton charge 
radius above 3 fm 2 where the contributions of the Dirac form factor and Foldy term to the 
neutron charge square radius tend to cancel. At the physical nucleon scale the value of r\ n 
is strongly dependent on the choice of quark spin coupling scheme, a consequence of the 
intrinsic relativistic nature of the constituent quark motion inside the nucleon. 
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FIG. 1: Light-cone time-ordered diagrams for the nucleon electromagnetic current. The gray blob 
represents the spin invariant for the coupled quark pair in the effective Lagrangian, Eq.@. The 
black circle in the fermion line represents the action of the current operator on the quark. Diagram 
(la) represents J^ N , Eq.(jHJ). Diagram (lb) represents J^ N , Ea. (|lU|) . Diagram (lc) represents J^~ N , 
Ea. ((TTj) . Diagram (Id) represents J^ N , Eq. ((T2)) . 
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FIG. 2: Neutron charge square radius as a function of the neutron magnetic moment. Results for 
the Gaussian wave function with a equal to 1 (solid line), 1/2 (dashed line) and (short-dashed 
line). Results for the power-law wave function with a equal to 1 (solid line with dots), 1/2 (dashed 
line with dots) and (short-dashed line with dots). Results for modified Gaussian wave function 
with a equal to 1 (full square), 1/2 (full diamond) and (full triangle). Results for modified power- 
law wave function with a equal to 1 (open square), 1/2 (open diamond) and (open triangle). 

1 14 
Experimental data from Ref.jg]. 
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FIG. 3: Neutron charge square radius as a function of the proton charge square radius. Theoretical 
results labeled as in figure 2. The experimental data points come from the measured value of the 



neutron charge square radius |gl, and from the experimental values of the proton charge square 
radius from Refs. 

, Q and Q, which are represented by the full circle, open inverse triangle 
and full inverse triangle, respectively. 
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FIG. 4: Individual contributions to the neutron charge radius related to the Dirac form factor 
Fi n (q 2 ) (r 2 n ) and from F2 n (q 2 ) (r| n ) as a function of the proton charge square radius. Theoretical 
results (squares, diamonds, triangles) for v\ n labeled as in figure 2. Theoretical results for the 
Foldy term, r| n , are bounded by the thick solid lines. 
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FIG. 5: Neutron charge square radius and r| n as a function of the proton charge square radius. 
Theoretical results for the Gaussian, modified Gaussian and modified power-law wave functions 
with a equal to 1, 1/2 and 0, labeled as in figure 2. Experimental points included in gray circle 
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